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Abstract 



We revisit the energy transfer necessary for the warped reheating scenario in a two-throat geometry. We study 
KK mode wavefunctions of the full two-throat system in the Randall-Sundrum (RS) approximation and find 
an interesting subtlety in the calculation of the KK mode tunnelling rate. While wavepacket tunnelling is 
suppressed unless the Standard Model throat is very long, wavefunctions of modes localized in different throats 
have a non-zero overlap and energy can be transferred between the throats by interactions between such KK 
modes. The corresponding decay rates are calculated and found to be faster than the tunnelling rates found in 
previously published works. However, it turns out that the imaginary parts of the mode frequencies, induced 
by the decay, slow the decay rates themselves down. The self-consistent decay rate turns out to be given by 
the plane wave tunnelling rate considered previously in the literature. We then discuss mechanisms that may 
enhance the energy transfer between the throats over the RS rates. In particular, we study models in which 
the warp factor changes in the UV region less abruptly than in the RS model, and find that it is easy to build 
phenomenological models in which the plane wave tunnelling rate, and hence the KK mode interaction rates, 
are enhanced compared to the standard RS setup. 
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1 Introduction 

String theory provides several promising and in many ways natural avenues to address the origin and properties 
of inflation as well as the origin of hierarchies of energy scales. For example, when the 6-dimensional internal 
space consists of several (perhaps many) long throats joined by bulk physical processes deep in a particular 
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throat will have their energies redshifted by an amount that depends on their position in the throat. In this 
manner one can obtain large hierarchies simply by cleverly choosing the internal geometry and positions of 
various physically important objects within it [Q: for example, placing the Standard Model on a stack of D 
branes deep inside a long throat can lead to a desired low SM energy scale [21 E]- 

An interesting class of inflationary models in such a multi-throat background is driven by the potential 
energy of D3 and anti-D3 branes J?. While the D3 branes are mobile, the anti-D3 branes are stuck at the 
bottom of one of the throats; during the "slow-roll" phase of inflation, the D3 branes move through the bulk 
and the appropriate throat (s) toward the anti-D3 branes OEI- Once they get within about a string length, 
a tachyonic mode appears and the branes and antibranes annihilate, ending inflation. The annihilation is 
expected to produce a large number of very massive closed strings j7j that will quickly decay into the lowest 
string states namely the (10-dimensional) graviton, the gauge fields and the dilaton, and their fermionic 
partners. 

Typically the annihilating branes are assumed to be located in a moderately short throat, so inflation 
happens at the scale required by CMB observations, about lO^'^-lO"'^^ GeV |9J. While the annihilation products 
have in principle enough energy to reheat our universe, the energy must first be transferred to our vicinity in 
the internal space, which usually means tunnelling through the bulk joining the throats. This issue received 
a good deal of attention recently [HI [HI ^1 , with the general conclusion being that generically the 

tunnelling seems to be too slow, but can be made fast enough for certain specific choices of the parameters in 
the model. 

Our aim is to make more precise several qualitative arguments and rough estimates used in the works cited 
above (the first attempt in this direction appeared in [Sj). Further, we investigate various mechanisms that 
may enhance (or, in some cases, suppress) the rate of energy transfer between the annihilation (A) and the the 
Standard Model (SM) throats. Following [5] [21 El; we will (mostly) work with two 5-dimensional AdS throats 
of constant (and equal) curvature separated by a UV brane (this setup can also be viewed as a two-fold copy of 
the Randall-Sundrum model [2], and we will occasionally refer to it as the doubled RS model). 

We begin in Section [21 by reviewing the setup of [HI El that will be the background for our tunnelling and 
decay calculations. In Section [3 we consider, in full detail, the wavefunctions of KK modes in this model and 
compare them to the wavefunction of the 4-dimensional graviton (zero mode). We find that the KK mode 
wavefunctions are enhanced only within narrow "resonance" bands; the width (in energy) of these bands can 
also be interpreted as the tunnelling rate from the A throat into the SM one. In Section IST^ we find that the 
narrow width is also responsible for the (perhaps surprising) fact that for tunnelling to actually happen, the 
SM throat must be much longer than the A throat. 

Even if KK modes localized in the A throat do not tunnel, their wavefunctions have a small tail in the SM 
throat. Thus they can decay into KK modes localized in the SM throat. In Section^lwe study the corresponding 
decay rates. Our results indicate that the decay into KK modes is faster than the tunnelling rates calculated 
in [SI [HI El and can dominate even if the SM throat is not as long as would be required for tunnelling. However, 
it turns out that the fast decay rates are not physical: in Section [SI we discuss modifications of the decay rate 
calculations due to different decay rates in the A and SM throats.^ We find that having different decay rates in 
the two throats changes the shape of the wavefunction dramatically: the mode will penetrate the side on which 
it decays faster only to a certain distance that is inversely proportional to the decay rate. This implies that 
long throats cannot enhance decay rates indefinitely; in fact, it turns out that the self-consistent decay rate is 
limited by the plane wave tunnelling rate considered in [3 [HI ■ We conclude that the plane wave tunnelling 
is realized, but only because the tunnelled particles decay in the SM throat fast enough. 

Lest the reader thinks that this whole exercise simply re-derived a simple result in a complicated way, the 
calculation of the self-consistent decay rate allows us to check that the tunnelling and decay are not slowed 
down further by the effect of complex frequencies (this possibility was pointed out in |l(Jj). 

In Section[Hlwe look at a possible mechanism for energy transfer that arises during SM throat relaxation after 
inflation. Intuitively it can be understood as follows: as the length of the SM throat increases, the spectrum 

^We were motivated by IIUI who argue that differing decay rates can lead to a significant suppression of tunnoUing. 
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of masses of its KK modes shifts down, while the A throat KK mode masses remain unchanged. When an SM 
throat eigenvalue approaches an A throat eigenvalue, instead of crossing they "repel" each other^: what was a 
KK mode localized in the A throat becomes a mode localized in the SM throat and vice- versa. Thus, KK modes 
and their energy could be "sucked" out of the A throat and into the SM throat where they could decay before 
another level repulsion would push them back. Unfortunately it turns out that the presence of an imaginary 
part of the mode mass (due to the decay) allows modes to avoid each other in the complex plane, and the mode 
switching between throats will not take place. 

In Section 13 we consider models that exhibit enhanced tunnelling rates, either due to resonant effects (such 
as the ones mentioned in |12 , 13 ) or because the shape of the potential barrier separating the throats is modified 
(and the barrier is lowered). The basic idea is to modify the throats near their top, where the effective potential 
barrier is largest. Examples we look at include two throats with different curvature radii, resonant tunnelling 
through the gravity box of jl5j . and a toy 5-dimensional model representing two throats joined smoothly in 
the bulk, in which the curvature is changing stepwise due to the presence of additional 3-branes.'^ We show in 
Section that low curvature in the central bulk region of the geometry can, under fairly general assumptions, 
lead to a drastic enhancement of tunnelling rates for the low-lying KK modes. We end in Section |51 with 
discussion and conclusions. 



2 Review: Model of a two-throat compact ificat ion 

We start by reviewing the 5-dimensional model of [5| I14| . We first summarize the background geometry, then 
concentrate on the fluctuations of the 5-dimensional metric as in |15| . 

2.1 Background 

We start with an infinite Standard Model (SM) throat and in later sections generalize our analysis to the case 
of a finite SM throat". Thus, our setup is Einstein's gravity in 5 dimensions parametrized by (a;^,z) with 
II = 0, ..,3, with a "Planck" (UV) 3-brane embedded at z = and an "annihilation" (A) brane embedded at 
z = —ZA < 0. The Planck and A branes have tensions V and —V {V > 0), respectively. The coordinate z thus 
runs from —za to oo. To simplify notation, however, we will use positive z values on both sides of the UV brane 
and will distinguish them by superscripts A (annihilation side) and S (Standard Model side). The action is 

J d-^xdzy/^{2MlR- A) - J d*xy/-g4{z = 0)V + J d^x^^g^{z = za)V. (2.1) 

kX, z — ZA we impose, for the mode analysis, Z2 boundary conditions on both the background and the fiuctua- 
tions. In Section [7| we discuss the appropriate boundary conditions for the tunnelling calculation. 
The background consists of AdS bulk with metric 

ds^ = cr2(z)(dx2 +dz2), (2.2) 

with the warp factor a(z) given by 



^This is known as level repulsion in standard Quantum Mechanics 

^Our term "3-branes" in this sense simply means domain walls with arbitrarily chosen tension; they are not supposed to represent 
string theory D3-branes. 

* According to a conjecture of 1111 . during inflation the SM throat warping cannot be stronger than the square of the A throat 
warping; it is assumed that after inflation ends the SM throat will settle into its "vacuum" (long) states, though this process may 
not be finished at the time relevant for KK mode tunnelling and decay. 
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and the curvature scale k (equivalently, the radius of curvature L = 1/fc) is determined by the cosmological 
constant A > as 

As is well-known 0, to get flat branes localized at definite points in the z direction, their tension V must be 
tuned to 

V = ±24M|fc , (2.5) 

where the plus and minus signs apply to the Planck and inflating branes, respectively. An extension of this 
formalism to the case where the curvatures on both sides of the brane are not the same is discussed in Section f7.1l 

2.2 Fluctuations 

We would like to study linearized fluctuations of the metric (|2.2() along the lines of 15 . We restrict ourselves 
to the metric fluctuations in the directions parallel to the branes and parametrize the full metric as 

ds^ = (r/^^ +\^)dx''dx'' + dz^ , (2.6) 

with the fluctuations = cr^/i^jy satisfying the 4-dimensional transverse-traceless gauge conditions. Setting 

h^.=e'P^a^/^ij^(z)e^,, (2.7) 

with — —p^ being the four-dimensional mass of the fluctuations, the linearized Einstein equations for the 
metric fluctuations reduce to a one-dimensional problem 



■^d'.+V{z) 



i'miz) = '^m'^ip„i{z), (2.8) 



with the potential V{z) given by 



The solution of H2.8() away from the branes is 

^^{z) - {mzfl'^ \Y2(Tn~z) + J2(mz)] , (2.10) 

where we have introduced 

l = z+\ (2.11) 

k 

to make the notation more compact. 

We now take ■0^ to be of the general form H2.10|l . with coeflicients TV;^, and iV^, on the SM and A 
sides of the Planck brane, respectively. The coefficient will be determined by the jump condition coming 
from (|2.8|l at the A brane, while and the ratio N^/N^ will be determined by the jump condition and 
continuity at the Planck brane. 

The following relation, derived using standard Bessel function identities and valid away from the branes, 
will be useful: 

d^-ipr,i{z) 3fc Fi(mz) -I- (3„iJi(to2) ^oioa 

^— — = a[Z) + m — -. (2.12) 

Tpm[Z) 2 Y2[mz) + Q,nJ2(mz) 

Equation (|2.8|l and the Z2 conditions imposed on the fluctuations hf^^ at the A brane lead to 

3fc 

dz^m{z Z^) = —(j{zA)'>Pm{zA), (2.13) 
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giving 



Qt = -^^y (2.14) 
Ji{mzA) 



Similarly, 1)2. 8|l at the Planck brane implies 



1 QjU 

2 (9.V^™(0+) - d,^^{0-)) = y Vm(0) (2.15) 



leading to 



g _ yi(m/fc) [Y2{mlk) + g4J2(m/fc)] + ^^(m/fc) [Y,{m/k) + Q^, Ji(m/fc)] 

J2(TO/fc) [Yi{m/k) + g4Ji(m/fc)] + Ji(TO/fc) [Y2{m/k) + g4J2(TO/fc)] ' ^ ' ^ 

Since the SM side is semi-infinite, we have a continuum of plane-wave normalizable modes with normalization 
constants given by JS] 




Imposing continuity of ipm at z = implies 

A _ .rs Y2{m/k) + Qg, J2(m/fc) 

'^Y2{m/k) + Qi.h{ra/ky ^ ' ^ 

2.3 Energy and curvature scales 

Before proceeding with our calculations, we should understand the physical scales in the problem. In the RS 
scenario that serves as our toy model, we have several mass and length scales: the 5-dimensional Planck scale 
Ms, an independent scale k that sets the curvature radius of the AdS space and the tension of the branes, the 
length ZA of the A throat and an (apriori independent) scale Ma of the brane physics that produces the KK 
modes which are, in the warped reheating scenario, responsible for the transfer of energy to the SM side. 

If our setup is to serve as a model of string theory compactifications, the parameters of our model must be 
matched to the appropriate scales in string compactifications. Hence we would like to take k to be somewhat 
below the strings scale Ms 

(2.19) 

19s 

with 7 a numerical factor of order say 10 that roughly corresponds to the compactification radius of the internal 
manifold, and gs the string coupling. Parametrized this way, the factor 7 is related |SlllU| to the flux quantum 
numbers of the Klcbanov-Tscytlin throat ^| E| as 

7-(Mi^)i/^ (2.20) 

The bare tension of the branes {D2> branes in string theory) is ~ M^/qs, with gs < 1. The brane scale Ma 
would be given by the tension of the annihilating branes times the warp factor at the position of the annihilating 
branes, 

Ms 
9 s 

9 s' 

« 7-, (2.21) 

ZA 
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where we have denoted a a = ^{za)- Thus, the effective mass scale on the A brane is somewhat (an order of 
magnitude, say) larger than 1/za- The 4- and 5-diniensional Planck scales M^^Mf, are related to the string 
mass Ms as 

Mlk = M| - Af3c,j3/4^5 , (2.22) 

where we have assumed that the compactification volume of the internal 5 dimensions is ~ . Expressed 
in terms of the string compactification parameters, the ratio kjM^ is 

The inflation scale iJ is different from Ma JIIj ^^^d in this setup is 

B ^ ^^M,a\ 

- -fk<T\ , (2.24) 

i.e., it differs from Ma by an extra factor of roughly^ cr^. 

Assuming H/M4 ~ (10~^ — 10^^) and the 4-dimensional Planck scale M4 being not very different^ from the 
string scale (and thus also from k) leads to the annihilation warp factor a a ~ (10^"^ — 1) To exhibit issues 
peculiar to warped inflation, in the remainder of this work we will assume that cr^ is in the low end of that 
range. 



3 Two throats full of modes 

In this section we study quantitatively the wavefunctions of KK modes likely to be produced in the decay of 
closed string remnants of the brane annihilation. As opposed to |S1 IIU) who considered modes in each throat 
separately, we study the modes of the full two-throat system. 

In Section [3. II we revisit the enhancement of the KK mode wavefunctions compared to the graviton deep in 
the A throat. Most of our conclusions agree with the ones that arrived at. The calculations are quite 

technical, so we relegate the details to Appendix IXI and provide only a summary in the main text. We find that 
most of the modes of the full two-throat system are suppressed at the A brane (compared to the graviton zero 
mode), even when their masses are of order of I/za] only modes that fall within narrow "resonance peaks" of 
width 

Smr^myk'^ (3.1) 

have their wavefunctions enhanced. These narrow peaks correspond to the discrete spectrum of modes of the 
A throat when the latter is viewed as a system on its own. 

The mode analysis allows us to exhibit tunnelling between the throats as a wavepacket decoherence effect in 
the combined system, which we describe in Section IX^ Wavepackets localized in the A throat will be composed 
predominantly of modes falling within the narrow resonance peak; hence the tunnelling rate T can be viewed as 
the rate at which the wavepacket decoheres and the destructive interference in the SM throat disappears. The 
rate is given precisely by the (energy) width of the wavepacket, 

m 

r^Sm^^. (3.2) 

This expression is analogous to the (inverse of the) "throat interaction time scale" given in JO] . 

In Section r3.3l we consider modifications of the tunnelling rate that arise when one makes the SM throat finite, 
with coordinate length zg. The spectrum will be quantized roughly in units of l/zs (for zs 3> za). We find 

^Derivation of 12.241 again assumes V5 ~ 1/k^. 

= Ms-f'^/gs', with standard assumptions M4 will be a few orders of magnitude larger than AIs 
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that for very long SM throats the resuhs are, as expected, unchanged. On the other hand, when the SM throat 
becomes so short that the mode spacing 1/zs becomes longer than the width Sm of the resonant peaks (|3.1|) . a 
mode localized in the A throat will remain localized there forever (at the linearized level; interactions can lead 
to decay of the mode). The surprising observation is that this suppression of tunnelling will occur when 

zs < ZA^ , (3.3) 

Interestingly, argue that during inflation, the SM throat cannot be longer than roughly ZAk/m which is 
much shorter than the above bound; after inflation ends, the SM throat presumable gradually relaxes to its full 
length, but the tunnelling may still be suppressed at the time the KK modes decay. 
We now turn to detailed calculations. 



3.1 Deep down in the throat: KK modes vs. the graviton 

In any model of inflation it is important that after inflation ends, the fraction of energy that ends up in 
gravitons be small. In brane inflation this means that production of gravitons from brane annihilation should 
be suppressed compared to other modes (in this case the KK modes) that can transfer the energy released by 
annihilation to Standard Model matter. 

The results of [H| indicate that the flat-space production rates of gravitons and KK modes in the decay of 
highly excited closed strings |Zj are roughly (up to a factor of order 1) the same. Hence, in the warped case, 
the relative production rates of KK modes and gravitons will be determined predominantly by the size of their 
respective wavefunctions. Thus, we should compare the wavefunctions of low-lying KK modes (m ^ k) to the 
wavefunction of the zero mode. The KK modes were analyzed in the previous section, while the zero mode, 
including its normalization, is 

\/k 

^0 = ■ (3.4) 

To quantify the relative magnitudes of the wavefunctions, we will study the ratio 

RA^) = —( — r> (3-5) 

and its physically meaningful integral (since the KK modes have a continuous spectrum) 

p{m)= / dm'Rlim'). (3.6) 

The function p{m) can be interpreted as the ratio of probability densities to find a KK mode of mass m' < m 
and to find a graviton at the A brane. If the couplings of such KK modes and the graviton to the closed string 
remnants of brane annihilation is similar, p{m) will roughly represent the relative abundance of KK modes vs. 
gravitons in the products of decay of the closed strings. Phenomenology requires that p{m) be large. 

One can consider KK modes with masses in two regimes: the first one will encompass "very light" KK 
modes in the sense that mzA <C 1- The second, "intermediate", regime contains masses such that mzA 1 
while mk <^ 1. In a full two-throat system we have a (semi-) continuum of modes in both regimes; on the other 
hand, an approximation with two separate, weakly coupled throats indicates that only a restricted subset of 
modes should be important in the A throat, namely the modes of the full system that approximate the modes 
of the A throat separated from the rest. This turns out to be the case. 

Very Hght KK modes, mzA ^ 1, are expected to be suppressed on the A side simply because the A throat 
is shorter than the SM throat (which we assume to be infinite here). The mass is so small that the volcano 
potential H2.9|l at the inflating brane is already large enough to suppress the mode: For kzA S> 1 the potential 
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is V{za) ^ '^^'^ would expect that for the KK masses m < 1/za the corresponding mode would be 

suppressed on the inflating brane. Explicit calculations in Appendix I A . 1 1 bear this expectation out: we find 

(3.7) 

and 

pM~2^«1- (3.8) 

Hence, production of such low-lying KK modes will be greatly suppressed compared to the production of 
gravitons. 

The situation is more interesting in the intermediate regime mzA ~ 1. In Appendix I A. 21 we find that for 
generic masses, the corresponding KK modes are still suppressed, 

2 

R^^ (3.9) 

so 

pirn) ^ — ^a\. (3.10) 

Here the equality holds only up to factors of order 1, which we denote by the ~ symbol. This suppression 
is counterintuitive: one would expect the Bessel functions entering the KK mode wavefunctions to begin ap- 
proaching their asymptotic regime and hence be enhanced (in the integrated ratio p{m)) by a factor of roughly 
<Ja^- To resolve this discrepancy, it helps to look at the plot of as a function of the mass m. Fig. E^b). The 
plot shows that modes with "generic" mass are, in fact, suppressed as indicated above; however, the plot also 
shows narrow resonant spikes in which the ratio i?^ is large. The width of the spikes is estimated in Appendix IXI 
and turns out to be given by the formula H3.1|l . 

(5m~mVfc^ (3.11) 

The contribution of the spikes to the integrated ratio p{m) turns out to be large and in fact confirms the intuitive 
argument for the KK mode enhancement: 

p{m) spike a-^^ ■ (3-12) 

This resonance phenomenon may seem surprising, but in fact it is to be expected in a system that is effectively 
composed of two weakly coupled potential wells. Those modes of the full system that have a small overlap with 
the "native" modes of the A throat (viewed as a separate system) will be suppressed there; the only modes of 
the full system that can be large are those that have a large overlap with the native modes of the A throat. 



3.2 Tunnelling as a decoherence effect 

Having modes of the full system allows us to exhibit the tunnelling of KK modes from the A throat to the SM 
one as a decoherence phenomenon: creation of a KK mode in the A throat is viewed as setting up a wave packet 
localized near the bottom of the A throat. As usual in quantum mechanics, to study the time evolution of the 
packet one can, e.g., expand it into eigenmodes of the full system (whose time evolution is known precisely). 
The initial expansion into eigenmodes is quite special {i.e., coherent), giving zero throughout most of the full 
system (namely in the SM throat). As the eigenmodes evolve with different frequencies, the evolution will 
destroy the coherence and the wavepacket's wavefunction will develop a nonzero generic amplitude in the SM 
throat, corresponding to the KK mode tunnelling out of the A throat. Let us call the decoherence time, i.e., the 
time it takes the (initially highly correlated) expansion coefficients to become essentially random, ttun ■ One can 
estimate ttun as the inverse of the spread in frequency, {dm)~^, of the initial wave packet, or the decoherence 
(tunnelling) rate Ttun as the spread itself. As we have argued in the previous subsection, the KK modes in the 
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A throat will be predominantly produced in the narrow bands around the peaks of ; the width of these peaks 
was estimated in (jA.23|) and we have 

4 1 4 

t^un = Tt™ ^ Jif^^ fcT™ ■ (3.13) 

The same expression was derived in Section 3 of [S], and a similar one was given Section 5 of |1L)| as the 
"interaction time scale"; however, the authors of ^U] assumed the UV scale to be the 4-dimensional Planck 
mass M4; our calculation above (and the calculations of 15) confirm that the relevant UV scale here is the 
curvature scale k which is expected to be somewhat smaller than M4. 

3.3 (Non-) Tunnelling into a finite SM throat 

In deriving H3.13|l we have assumed that the SM throat is infinite, even though in reality it will be long but 
finite; we will now clarify how long it must be to be for our preceding calculations to apply (inasmuch as our 
results agree with the ones in this restriction applies to those works as well). Hence, in this section we 

will assume the SM throat is cut off a,t z = zs > za, with the > sign meaning at least 1 order of magnitude. 

If the SM throat is finite with length zg, KK modes of the full system will be quantized roughly in units of 
SsTn l/zs- One might think that simply having zs > za would mean that from the A throat point of view we 
a have a quasi-continuous spectrum, but this is not the case. The fastest way to see it is to consider the resonant 
spikes in discussed in Section im and Appendix I A. 21 the width of these spikes is (5m ^ {m/k)'^z^^ <C z^^. 
In fact, if the spacing Ssin is larger than the width of the spikes, it is likely that none of the modes present in 
the system will fall within the interval Sm, indicating a suppression of tunnelling. This argument is a bit naive 
in that in addition to the J^m spaced modes (that one can consider to be native to the SM throat) there will 
be extra modes with spacing SAfn ^ ^/za that are mostly localized in the A throat; however, precisely because 
these modes are localized in the A throat, wavepackets consisting mainly of such modes cannot tunnel. 

To put some firm math behind the intuition we consider modes in the finite two-throat system explicitly. 
The boundary condition H2.14|l now applies to the SM side as well (with the obvious modification), 

Qi-~^^y (3.14) 

Ji(mzs) 

Modes of the full system have masses m for which and satisfy H2.16|l . 

The discrete spectrum of modes implies a change in normalization: just like the zero mode, the massive 
modes are now normalizable in the conventional sense, 

dz^Pf^{z) = l. (3.15) 

Fortunately, the normalization integrals for the wavefunctions given by H2.1()(l can be evaluated in a closed form. 
We list them here for future use and define a shorthand notation: 

Cj{x) = J dxxJlix) - [Jl{x) - Ji{x)J^[x)\ , (3.16) 

Cy{x) = j dxxYiix) = [Yiix) - Yi{x)Y3{x)] , (3.17) 

Cjy(x) = J dxxY2{x)J2{x) = [2Mx)Y2{x) - Mx)Y3{x) - Y,ix)J3ix)] . (3.18) 

To get oriented, let us first look at modes that are mostly localized in the SM throat (these are the equivalents 
of the "generic" modes discussed in Appendix IA.2|I . Such modes have ~ 1 {i.e., the Bessel functions 
entering H2.14|l have generic values); this implies 

Q™^-^»1, (3.19) 
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meaning that the masses of such modes he near the zeros of Ji{mzs), confirming the intuition that the mass 
would be quantized in units of 1/2:5. Since 3> Q^, the normahzation wiU be determined predominantly by 
Cj{mzs) and works out to be 



TO m? 



Ng ~ , , ~ -pf— . (3.20) 

Q^Cj{mzs) fc2^ 



Continuity of iprn, Eq- H2.18|) . then gives 
and for the wavefunction itself we find 



Comparing this with the zero mode wave function (|3.4|l we find 



N^^Ng, (3.21) 



^^{za)-^:^. (3.22) 



^^^"^^ 1^ (^-^^^ 
Summing over all modes whose m ^ za we find the equivalent of IIA.12|l , 

p{m)^ J2 i?|(m')~^. (3.24) 

m' <m 

(We have not explicitly considered the very light modes, but just as in the continuous case, their contribution 
scales the same with m/k.) This means that the generic (or SM throat) modes are small in the A throat. In 
the continuum case we found that most of the KK mode presence was due to a narrow range of modes that we 
called resonant. In the discrete case such modes occur when is (nearly) zero and large is large: = 
implies, via H2.16|l . 

Qt^^^l. (3.25) 

In this case the biggest contribution to the normalization integrals can come either from Cy{zs) ^ mzg or from 
Cj{za) ~ mzAk'^/m*, depending on the relative size of za, zs and 1/fc. 
If the SM side dominates, namely for 

ZS > ZA^ , (3.26) 

TO^ 

the mode will be mostly localized on the SM side despite being large in the A throat, simply because the SM 
throat is so long that the probability of finding the particle in it is close to 1. This is the case that closely 
approximates the continuum setup: the resonant KK modes are large in the A throat, yet if a wave packet 
localized in the A throat is formed from them, it will tunnel to the SM throat as discussed in Section 13.21 
Eq. H3.26(l confirms the intuitive arguments given at the start of this section: precisely when the former is 
satisfied, the spacing of the SM throat modes is fine enough that some will fit into each resonant spike. 
On the other hand, if the A side dominates, namely for 

k^ 

ZS < ZA— , (3.27) 

TO 

we find a surprising result: the KK modes that are large in the A throat have a very small probability of being 
detected in the SM throat, and if a wave packet of such modes is formed in the A throat, it will generically 
remain there^; only a small part of order zsm'^ / {zAk'^) < 1 of these KK modes will tunnel. Note that this can 

^It may be possible to form a coherent wave packet that will appear very briefly in the SM throat, but it will decohere quickly 
and the particle{s) will reappear in the A throat again. 
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be the case even if the A throat is much shorter: £.17., if m ~ \Iza ^ 10~'^fc, the A throat can be ahnost 12 
orders of magnitude shorter! 

The preceding arguments provide a further constraint on the relative length of the throats: since the modes 
relevant to tunnelling have masses of order 1/za, Eq- H3.26|l implies 

zs > a\zA . (3.28) 

Written in terms of warp factors, we have a constraint on the Standard Model warp factor as = '^{zs) that 
reads 

as<a\. (3.29) 

This is a stringent constraint; in fact, |llj suggests that during and shortly after inflation, the warp factor of 
the SM throat is no less than the square of the A warp factor, 

<ys'^'JA, (3.30) 

and hence the inequality (|3.29() would not be satisfied.® 

In Appendix ^ we show that a condition similar to H3.26|l can be derived in a toy model of a (5-function 
potential barrier separating two finite flat boxes. A natural question to ask is what is the physical interpretation 
of the factor /k^ relating the two throat lengths in H3.26|l and its analog in the (5-barrier case. We propose 
that, in general, tunnelling will occur only if the level spacing of the target side is finer than the level spacing 
on the source side multiplied by what one may call "plane wave tunnelling probability" P, 

{5m)target < P{Sm) 

source ■ 

(3.31) 

The tunnelling probability P is the "usual textbook" probability of a particle that starts as a plane wave 
far from the barrier on the source side to emerge as the corresponding plane wave on the target side. For 
the Randall-Sundrum model P was calculated by and found to be P ^ m*/k^; we show in Appendix IbI 
that (|3.31l) also applies to the S potential barrier. 

In retrospect, at least naively, this looks like an obvious statement: the size of the wavefunction squared, i.e., 
the probability density, in the target well is suppressed by the tunnelling rate with respect to the probability 
density in the source well; hence their lengths must differ by the corresponding factor for the probabilities (which 
are integrals of the densities) to be comparable. However, the statement uses the tunnelling rate of plane waves 
(which is independent of the lengths of the wells) to compare the tunnelling probabilities of wave packets made 
of standing waves in finite wells (and in finite wells the spectrum and tunnelling probabilities do depend on the 
lengths of the wells), so perhaps the result is not so obvious after all. 

Incidentally, the above result, if true in general, would be applicable also to the multibrane setups considered 
in Section[71 the enhancement of the tunnelling rates of plane waves implies that the ratio of throat warp factors 
need not be as large as indicated in (|3.29|) . 

Before we end this section, a comment is in order about an apparent discrepancy with 0, whose authors 
claim that tunnelling rates are essentially {i.e., up to factors of order 1) unaffected by making the SM throat 
short. The calculation in Appendix A of [S] appears to agree with ours above (although it is formulated in 
different terms), but their interpretation is incorrect in the sense that pi split standing waves into in- and 
out- going parts and then regard the latter as independent, interpreting the outgoing component as representing 
a flux of tunnelled particles. 

4 Energy transfer via decay without tunnelling 

If the throat lengths are such that A throat KK modes remain localized in the A throat indefinitely, the said 
KK modes will decay via the nonlinear interactions present in the Einstein-Hilbert Lagrangian (and via the 

*As wo explain in Section 17.21 it is possible, if somewhat contrived, to relax I3.29i if the throats have different curvatures. 
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coupling of the 5-dimensional gravitational fluctuations to the SM brane fields). One might expect that for the 
lowest A throat KK mode the only decay channel would be pairwise annihilation into gravitons, but this is not 
so: the wavefunction of A throat KK modes has a small tail in the SM throat, and the tail actually causes a 
sizable decay rate into SM throat KK modes. In this section we study this decay rate and find that in general 
it appears to be faster than the plane wave tunnelling rate. 

Before we dive into the calculations, however, a cautionary note is in oder. It is clear that the amplitude of 
the A throat KK wavefunction tail in the SM throat depends to a large degree on the potential shape in the 
UV region. For example, Ref. fTS" found that tunnelling out of a full AdSs throat into the bulk is suppressed 
compared to the analogous quantity in the RS model; this suppression can be viewed as due to a modification 
of the effective potential entering the Schrodinger Equation \2.8\\ . While the Randall-Sundrum geometry is a 
good approximation deep in the throat, the region near the UV brane will presumably look very different in 
(a dimensional reduction of) a full string compactification. With this possible loophole in mind, let us proceed 
with the actual calculation of the decay rates. 

We start by quantifying the interaction strengths, following 0. The 5-dimensional Einstein-Hilbert La- 
grangian M^^/gR gives rise to an infinite series of interactions of the (canonically normalized) fluctuations h 

of the metric g around its background value g (g ~ g + M^^^h). Schematically they are of the form 



1 



M3n/2 — 3 
5 



V?V^ {g-y^'h^^ (4.1) 



where the derivatives can act both on the background metric and the fluctuations. We will only be interested 
in the powers of the warp factor and will neglect all order 1 coefficients (in particular, we will not attempt to 
keep track of the factors arising from various possibilities of contracting the indices in (|4.1(l ^. 

A z-derivative acting on g~^ will give kag^^ (derivatives in the directions are zero). A z-derivative acting 
on h will produce several terms that can be roughly characterized as being either of the form kah or mh, where 
m is the mass of the KK mode. Derivatives of h in the directions will simply multiply the wave function by 
the corresponding momentum. Hence the interaction (|4.1() leads to the following 3 basic types of terms, 

^ ^ ^{g~'Y^^hJ'{{p-p,mm),mka,{kaf] (4.2) 



(We consider the p ■ p, mm terms to be of the same type; note that the p ■ p and mm symbols should not be 
taken at face value; they could mean, e.g., pi ■ p2, mim2 if the derivatives act on different KK modes.) The 
terms differ by their z dependence and dependence on momenta. 

To obtain an effective 4-dimensional interpretation, we perform the integral over z in H4.2|l : the results will 
be an effective 4-dimensional Lagrangian. We will treat all terms in (|4.2|l together; denote the power of ka by 
a (so a = 0, 1, 2); then we can summarize (|4.2f) as 

3^/2-3 / d'^xdzy/g{g-y^^ hi...hnWi^a)k°'a°' (4.3) 

where u'(o) is a combination of mm and p ■ p, W(^ij = km and W(^2) = 

We recall from Section that the z-dependence of the graviton fluctuations h is given by the KK mode 
wave function and the appropriate prefactor (we neglect the Lorentz tensor structure), 

h{x, z) = a^^^{z)^{z)(l3{x) (4.4) 

where (pix) is the 4-dimensional wavefunction that is usually taken to be simply exTp{ipx). The n-point couplings 
can then be written concisely as 

£„ = I d'^xC,ri4>i{x)(j)2{x)4>:<,{x) , (4.5) 
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where the "couphng constant" C„ (which in fact can contain factors of momenta) is the overlap integral along 
the z coordinate, 

STTTiTl \dz^9 a"/2u;(„)fcV"Vi(^) • • ■ V'n(^) • (4.6) 

Amplitudes involving at least one graviton deserve a special consideration.^ From the perspective of the effective 
4-dimensional theory, (4-dimensional) general covariance requires that the graviton couple to the (4-dimcnsional) 
stress tensor, which must in this case be proportional to (4-dimensional) momenta; hence amplitudes that include 
a graviton can only be nonzero for a = 0. 

Evaluating the integral (|4.6|l is in principle (numerically) straightforward, but giving meaningful estimates 
is somewhat tedious because the dominant regions in the z-integrals depend on how many A throat KK modes, 
SM throat KK modes and gravitons enter the amplitude. Hence we will calculate the 3-point coupling (|4.6() for 
1 A throat KK mode with mass toq decaying into two SM throat KK modes with masses toi,2- 

In principle one could contemplate the coupling of the A throat KK mode to a graviton and an SM throat 
KK mode, or to two gravitons, but these couplings vanish: because of the graviton profile, a 3-point amplitude 
with a = and one graviton reduces to the scalar product (of the z-dependent wavefunctions) of the two 
remaining modes in the amplitude, which is zero whenever the two modes are different. (As noted earlier, the 
a 7^ couplings involving gravitons must vanish by 4-dimensional general covariance.) 

To perform the integration, we approximate the Bessel functions by their small argument approximations and 
asymptotic forms respectively, where appropriate. In the region of z < 1/ \/mk, which we will call the UV regime 
because it describes the region near the UV brane, the KK mode wavefunctions have the same functional form 
as the graviton wavefunction, but their normalization constants differ. For integrals that are mainly localized 
in the UV region we can neglect the oscillatory nature of the Bessel functions at large arguments; for integrals 
that receive large contributions from the asymptotic region, we will have to take the oscillations into account. 

Near the UV brane in the A throat (z < 1/^/mk), the wavefunction of an A throat KK mode behaves 
roughly as 

V;(z)^y^ayV/2(z), (4.7) 

while far from the UV brane {mz ^ 1) it is 

V'(^) ~ — = = Vka]/'^ = const . (4.8) 

\/ZA 



In the asymptotic region of the SM throat, the wavefunction of an A throat KK mode is 

V.(z)^^^ = Vfcay^^= const. (4.9) 
The wavefunction of an SM throat KK mode with mass in the asymptotic region of the SM throat is 

tljiz) ^ = Vka]/^ = const , (4.10) 

Vzs 



while in the A throat it is (here we assume niiZA < 1, so only the UV regime is realized) 



1 



- 1^^^^ = ■ (4-11) 

As we remarked above, the functional form of the KK modes in the UV regime and the graviton (given by H3.4|l ) 
is the same, only the normalization constants differ. 



^We thank Henry Tye, Xingang Chen and Irina Mocioiu for corrections, discussions and suggestions on this point. 
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Under the approximations H4.7II - H4.11|I the integrand in H4.3II wiU be a certain power of the warp factor a{z) 
(multiphed by various constant factors). If this power is negative or zero throughout most of the throat in 
which the integral is performed, the IR region of the throat wiU dominate the integral; however, if the power of 
(t{z) is positive, the UV region will dominate. 

The amplitude of interest here is the decay of an A throat KK mode. We will assume that we are dealing 
with the lowest-lying A throat mode; higher modes either decay into the lowest one or decay via the same 
channels as the one we are considering. The coupling to two SM throat KK modes is dominated by the IR 
region of the SM throat and is largest for the derivative coupling {a — Q). The result is 

yS c k^'"^ ™0 1/2 -3/2 I. 

In this case, the oscillatory nature of the wavefunctions cannot be neglected since the integrand is large in the 
asymptotic region. To remind us of the fact that the naive expression is an overestimate, we have inserted a 
factor of ^ in 14.12|l to account for the extra suppression. The factor ^ can be roughly estimated as follows. 
The integrand is ~ a~^^'^{z) times a factor of type cos(miz) for each wavefunction. The naive estimate of the 
integral, integrating z from to zg, would be ^ a^^^^/fc. An improved estimate can be made by noting that 
W cos(miz) can be written as a sum of sines and cosines of various sums and differences of the arguments rriiZ. 
Hence a better approximation to the integrand would be ct~'^/^(z) cos(mz); for a generic value of m the leading 
term in the integral would go as / dza~^/'^{z) cos(mz) ^ a^^"^ jm. Comparing this result with the naive value 
fixes the correction factor to be ^ ^ Xjimzs) asjoA- We note that if the masses entering the amplitude 
are tuned such that (say) toq = ^\ +^2, the naive estimate is actually correct (as can be seen by taking a limit 
m — > of the improved integral). However, this "resonance" cannot be used to enhance the coupling, because it 
happens precisely when the decay products are on the threshold and have zero available phase space. Including 
the correction factor gives an improved estimate of the 2S coupling, 

/-S ^^^"^ ™0 -1/2 -1/2 

a,2S - ^^o-^cT^ ' <Js- (4-13) 

The couplings can be written in various forms using mo/k ~ a a and defining the Standard model mass scale 
MsM = Mc,as', the ratio k/M^ (an independent parameter in the RS model) is in a string compactification 
determined by the compactification moduli. 

To obtain decay rates from the above couplings, one must perform a sum over accessible final states. In 
Appendix iniwe perform the sums in detail; however, up to numerical factors, the correct results can be obtained 
by simply assuming that all momentum factors are of order mg and simply multiplying the resulting coupling 
(or decay rate) by the number of accessible modes, which in case of a decay into two KK modes it will be n^j^. 
For the decay into two SM throat KK modes, the coupling H4.12|l and (jC.lOp give the decay rate 

TA,2S-k(^^^\\af. (4.14) 

This rate is quite fast; in particular, by making the SM throat long enough, it could be made much faster than 
the plane wave tunnelling rate H3.13|l . However, in the next section we show that the actual decay rate is slower 
than (jnijl . 



5 Tunnelling is back: Effects of decay on KK mode wavefunctions 

In the previous section we have found the decay rate H4.14|l of an A throat KK mode in the SM throat. The A 
throat KK modes can also decay in the A throat itself by pairwise annihilation into gravitons; the rate of this 
process was estimated in |2] to be 

rA,2g-ka\. (5.1) 
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Generically, the throat lengths are expected to be such that Ta,2S 3> TA,2g- In this section we would like to 
incorporate the effects of unequal decay rates into our treatment of the two-throat system. 

Let us first use intuitive arguments to show what we expect to find. Wavefunction decay can be modelled 
by the presence of a sink (an imaginary term in the potential). If the sink varies as a function of the spatial 
coordinate (in our case, z), the wavefunction will tend to decay at different rates in different regions. Such 
a process would increase its gradient energy too much; to counteract it, there will be a non-zero probability 
current from regions with a small sink to regions with a large sink. In the case of two throats with different (but 
constant within each throat) sinks we may expect that the throat with the smaller sink will act as a "reservoir" 
from which the probability density will slowly "seep" into the throat with the big sink. The seepage between 
throats is expected to be slow because of the presence of the potential barrier. In the throat with the big sink, 
we have a source (incoming current density) at the top of the throat while the sink is present throughout, so we 
would expect the wavefunction to be largest near the source and decrease as a function of the coordinate z along 
the throat. A toy model of such a situation for a 5-function barrier separating two square wells is presented in 
Appendix^ the model confirms the intuition outlined above. 

An additional feature of the RS setup is that the decay rates, and hence also the sinks, depend on wavefunc- 
tion overlaps; as we argued above, the wavefunction of the decaying mode will, in turn, depend on the decay 
rates. Thus a would-be large decay rate in the SM throat would mean that the overlap integrals are actually 
small, so the actual decay rate cannot be too fast. How fast is too fast? Physically, the decay rate in the SM 
throat cannot be faster than the rate at which probability density leaks from the A throat into the SM throat 
- the presence of the barrier here turns out to be the limiting factor in the decay rate. We show that in the 
presence of a large sink, the wavefunction in the SM throat actually becomes a (spatially decaying) outgoing 
wave, and the "leakage rate" becomes simply the plane wave tunnelling rate through the barrier. 

We now describe detailed calculations supporting the intuitive arguments just given. We start by analyzing 
the wavefunctions of the KK modes in the presence of sinks and then give an improved calculation of the decay 
rate Ta,2S from Sectional 



5.1 KK modes in the presence of sinks 



The effect of decay of KK modes will be represented by an imaginary term iS in the potential V{z) in the 
effective Schrodinger equation (|2.8(l . We will take S to be constant in each throat, but the values and 
in the A and SM throats, respectively, will be different (this is a simplification of the real problem where the 
decay term varies smoothly with z; the advantage is that it has a simple analytic solution). The 5-dimensional 
graviton fluctuations still have the form H2.7|l . but now we allow the frequency p° = = 771 + to be complex. 
The effective Schrodinger equation is then 



ldf + Viz) + iS 



i^m{z) = -Uj'^1pm{z), 



'4'miz) 



The solution is a generalization of 12.10|) : in the A and SM throats we have, respectively, 

N^,VkS [Y2{X'^S) + QiJ2{X'^S)] for z < , 

Ng^/kS [Y2{X^S) + QiJ2{X'^z)] for z > , 
where the (complex) wavenumbers A"^ , A"^ are now different from the frequency (and from each other) , 

A^ 
A^ 



2iSS. 



(5.2) 



(5.3) 



(5.4) 
(5.5) 



^"We were motivated by Kofman and Yi )10l who suggest, by treating the two-throat system as two separate throats coupled 
weakly through the barrier, that differing decay rates (represented as differing imaginary parts in the modes' frequencies) will 
suppress tunnelling rates. 
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The matching conditions at the ends of the throats and at the UV brane are direct generahzations of the ones 
derived in Section at the ends we have 

(5 6) 



(5.7) 



MXSSs) ' 

while the jump and continuity conditions at the UV brane can be written as 

A Y,{X^/k) + QMX^/k) ^ s YiiXVk) + QiMXyk) 
Y2{X^lk) + QiUX^Ik) Y2{Xyk) + Q'lMX'yk)' 

N:^,[Y2{X^/k) + QiMX^/k)] = Ng[Y2{X''/k) + QiMX''/k)] . (5.9) 

Let us see how the spectrum wih look like. We will assume (based on the decay rate calculations) that the 
magnitude of the sink terms is much smaller than the curvature scale k^ (indeed if they were not, bulk effects 
would invalidate our effective field theory analysis), though the sink terms can be comparable in magnitude 
to the typical (squared) KK mode mass scales l/z^. Then we have jA'^'"^! <C k and we can apply the small- 
argument approximations (jA.ip to the Bessel functions in the UV brane jump condition H5.8|l . Analogously 
to the sink-free case, the jump condition can only be satisfied when one of the Qm coefficients is large. For 
definiteness let us look at the case when becomes large (such modes will again be the A throat modes). 
From H5.6I) we find that Ji{X^za) must be small, so A'^i^ must be near one of the roots of Ji, all of which are 
real. Hence A^ will be (nearly) real (it will have a small imaginary part that is not important for us). This 
fixes As to be 

Xs^^J{X^r-2z{SS~SA), (5.10) 

so A5 = a + i6 will have a substantial imaginary part b. The most important consequence of b is that in the 
asymptotic region of the SM throat the wavefunction Y2{X^z) + Q^J2{X^ z) drops off roughly as as exp(— 6z). 
This can be seen from the asymptotic form of the Bessel functions, namely 

^pKxj2Jy(x) w cos(x- iy7r/2-7r/4), (5.11) 

^J'KxlTYy{x) w sin(a; - iy7r/2 - 7r/4) . (5.12) 

Indeed, taking into account H5.7(l we find 

2 

V^(^) ^ ^f|^cos[A«(5-%)] 

7rCOS(A-'2:5) 

- Nge-"'. (5.13) 

Note that remains finite in the limit Z5 — > 00. 

This is an interesting result: it says that as long as the wavenumber A;^ has an imaginary part (which will be 
the case whenever ^ S^), the A throat KK modes will effectively penetrate the SM throat only to distance 
at most Zp ^ 1/6, no matter how long the throat is. 

An important consequence of the finite penetration length, if it is shorter than the throat length zg, is that 
it can lead to a substantial reduction of the overlap integrals (I4.()|l that are the coupling constants in the decay 
rates. We will now argue that the decay rate H4.14|l is so fast that it implies a penetration length Zp much 
smaller than the length of the SM throat; hence, in a sense, the rate Ta,2S is inconsistent and we must find a 
way of estimating the physically realized decay rate. 

To obtain the penetration length Zp implied by the decay rate r^_2S, we need to find the sink that 
corresponds to This problem is solved in Appendix IeI the result is 



k^a\a^^ , (5.14) 



16 



so 



Clearly, unless k/M^ is extremely small, we find 



Zp<zs, (5.16) 



meaning, as we claimed above, that the penetration length is much shorter than the throat length and hence 
the physically realizable decay rate rA,2S will be smaller than the result H4.14(l . We now estimate rA,2S- 

5.2 Estimating the effective decay rate 

We have an A throat mode that decays by interactions that are localized in the S throat. We argued above that 
its wavefunction will decay exponentially as a function of (coordinate) distance in the S-throat. The situation 
can be described as a reservoir of probability density in the A throat that slowly leaks into the S throat, where 
the probability density decays. Clearly, the probability decay rate cannot be faster than the rate of leakage; in 
fact, since the probability does not accumulate in the throat either, the decay and leakage rates must equal. 
The leakage rate is given by the probability current 

jz^-rd^ (5.17) 
m 

evaluated at the top of the throat. The current (|5.17() will consist of two terms j = ji + j2, with ji proportional 
to the imaginary part of the coefficient Q^, while j2 is proportional to the imaginary part of the wavenumber 
A'^. Let us first look at ji: using the asymptotic forms of the Bessel functions and bzs > 1, we find 

Qi - (5.18) 

so the wavefunction in the SM throat is actually the outgoing Hankel function _ff^(Az)! The parameter Xz is 
complex, so the amplitude of the wavefunction decays exponentially, as we have shown in the previous section; 
near the top of the throat we find 

ji ~ ka% ma\ , (5.19) 

namely the rate at which plane waves tunnel between two infinite throats, considered in |14l IHl ll()| . The part 
j2 of the current can be estimated as 

Im 



J2 Ji 



Re A^^ 



(5.20) 



The negative sign is always present and means that the decay of the wavefunction further slows down the 
"leaking" of the probability density into the SM throat. An important question now is whether j2 is comparable 
to ji, that is whether the imaginary and real parts of the SM throat wavenumber A'^ are comparable. If it 
were so, the decay rate could be significantly suppressed even compared to the plane wave tunnelling rate (this 
possibility was also pointed out in \W\). As we argued earlier, the imaginary part of the wavenumber determines 
the penetration length Zp as ~ Im {A"^}. In the following subsection we will estimate the latter and show 
that, at least in our approximation, the current j2 is indeed negligible, so the effective decay rate roughly equals 
the plane wave tunnelling rate, 

^A,2S ~ ^tun ~ ma\ . (5.21) 

Thus, while the naive decay rate H4.14|l is not the physical one because it is larger than the tunnelling rate, it 
does have a physical meaning: it is in turn small enough that current j2 remains negligible compared to ji . 
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5.3 Estimating the penetration length 



We would like to calculate the penetration length Zp = b^^ of the A throat KK mode wavefunction into the 
SM throat self-consistently, at least within our approximations. The strategy is to calculate the decay rate 
Ta,2s{zp) as a function of Zp, and equate it to the maximum physical rate (|5.21|l . 

The calculation will follow the same line as in Sectional with the following modifications. First, we will 
assume that the decaying wavefunction, instead of decreasing exponentially with the characteristic length Zp, 
equals the naive wavefunction for z < Zp, and is zero for z > Zp] this means the overlap integral H4.6|) will only 
extend to z — Zp. The wavefunctions entering the integral are still approximated by H4.9|l and H4.10|l and the 
resuh is (cf. gj^) 

Ci2S - |-S^-o^ay^cTsa(z,)-^/^ . (5.22) 

5 

The correction factor ^ is now, in analogy with the discussion below (|4.12|1 . roughly 

^"~(T^V(zp), (5.23) 

SO the improved coupling can be written as 

L.l/2 

Ci,2S ~ -^woCTT<^S<yizp)-'^^ . (5.24) 

Calculating the decay rate also involves summing over accessible final states. Since the wavefunction of the 
decaying mode only penetrates the SM throat to distance Zp, it will only couple (appreciably) to SM throat 
modes whose mass is large enough for them to not be suppressed for z < Zp, so we will only include modes of 
mass m > 1/zp in the sum over final states. Roughly speaking this means that the number of accessible final 
states changes from zsniQ to Zpmo. Applying the results of Appendix ICl we find 

rA,2sizp) ~ ^^Ws'^-Hzp) ■ (5.25) 



Equating TA.2s{zp) with the actual rate (|5.2H) then gives 

^ - 1^4^ (5-26) 
a A [Ml all 

< 1. (5.27) 

Put in words, we find that the penetration length Zp is longer than the A throat length za- Hence the imaginary 
part of the wavenumber in the SM throat, Im {A'^} ^ Zp^, is much smaller than the real part Re {A"^} ~ z^^, 
which is what we wanted to show. 



6 Level repulsion and crossing during SM throat relaxation 

Recall that according to the conjecture of JI], the SM throat will be shortened during inflation by the effects of 
interactions of the A and the SM throat length modulus. After inflation ends, the SM throat should relax to its 
"full" length. While the quantitative details of the relaxation, such as its rate, are unknown, one can imagine 
modelling it via a time-dependent SM throat length zs in the doubled RS model that we employ to represent 
the full string theory setup. 

A changing SM throat length leads to an additional possible mechanism of energy transfer. As the length of 
the SM throat increases, it will periodically go through an interval where the lengths of the A and SM throats 
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are tuned in the sense that there is a (near) degeneracy between pairs of one level from each throat. As is 
usual in quantum mechanics, the energy levels of coupled systems avoid crossing each other, which in this case 
means that a level that started out as an A throat KK mode will, after passing through the tuning point, 
become an SM throat KK mode and vice versa. This mechanism could, if the SM throat relaxation is slow 
enough, "suck" the energy out of the A throat and into the SM throat, where the KK modes would quickly 
decay into lower- lying states that cannot switch back to the A throat. We start by presenting a simple version 
of this mechanism that neglects the presence of complex parts in the frequency and wavenumbers induced by 
mode decay; then we investigate whether properly accounting for the said complex parts changes the picture 
qualitatively. 

We will assume that the SM throat relaxation is adiabatic, at least as far as the A throat KK modes (and 
SM throat modes of comparable masses) are concerned. The usual condition for adiabaticity, w/w^ ^ 1 implies, 
via zguj ~ const, a condition i <ti crj^/ag. We will simply assume this condition is satisfied. Let us now look 
at how the KK modes evolve with a changing zs while = const. To get oriented, we first take the sink-less 
case of real A'^ = A'^ = m. The spectrum-generating matching condition H2.15|l can be re-written in a more 
symmetric form as 

Fi + QjJ, Fi + QjJi _ 
where the Bessel functions are all evaluated at m/k and the coefficients Q are given by 

Q^'--^f^y (6.2) 

Let us start with zs such that the two throats are detuned (that is, only one of the coefficients Q is large). 
Recall that an A throat KK mode has a large and a small (^ 1) Qfj^. As zs increases, the mode's mass 
remains approximately constant until Q^^ starts becoming large as well. At this point the mass m starts moving 
such that mzs will remain approximately constant, will become small (~ 1) and will become large - 
what started as an A throat KK mode has become an SM throat mode. The switch-over will repeat itself in 
reverse when zs decreases further and the throats become tuned again. Of course, if we had started by tracking 
an SM throat mode, it would switch over into the A throat at the first point of tuning, and then back to the 
SM throat at the second one. This is an example of a standard quantum-mechanical level repulsion: at each 
point where the two throats have tuned lengths, the levels of one A throat mode and one SM throat mode avoid 
crossing by switching their "home throat" . In principle, this switching can facilitate energy transfer between 
the throats; if it is slow enough to allow the modes to fully relax, a KK mode particle originally located in the 
A throat will relocate to the SM throat where it can decay to lower KK modes very efficiently. If the A throat 
mode is the low-lying one, it is quite likely that the products of its decay will have masses below the lowest 
lying A throat KK mode and hence will not be subject to switching anymore. 

This simple picture has an obvious loophole: when one takes into account the presence of the sink(s), the 
mode wavenumbers A"*'"^ acquire (different) imaginary parts and it becomes possible for the levels to avoid each 
other in the complex plane even though the real parts cross. It is easy to see that this will happen when the 
imaginary part of one of the wavenumbers is large: in that case the corresponding coefficient Qm cannot become 
large (as we show below), meaning that the mode can never switch to that side. On the other hand, one would 
expect that if the sink(s), and hence the imaginary parts of the wavenumbers, are very small, the system should 
behave similarly to the one without sinks at all (in other words, the equations governing the eigenvalues are 
analytic and we do not expect any discontinuities for S* — > 0). Hence there should be a maximum sink that still 
allows the mode switching to take place; above that sink the real parts of the wavenumbers can cross and modes 
remain in their respective throats. Let us find the condition for the switching to take place, and see whether 
our model satisfies it. 

The spectrum is determined by (|5.8|) . It is clear that as long as both the (real) mass m and the sink S are 
much smaller than the curvature scale k, the only way to satisfy H5.8|l is to have at least one of the coefficients 
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Qm large. Suppose we start with detuned throats and pick an A throat mode, i.e., is large: to leading 
order, 

^™ ^ X^J^Yi ^ A2 ' ^^-^^ 

where we have introduced a shortened notation for the Bessel functions by denoting Yf = Yi{X^ /k) etc, and in 
the last expression we have neglected the difference between the wavenumbers. As we have noted before, H5.6|l 
implies that A"^ must be (nearly) real. 

For the mode to switch throats, must become of the same order at the tuning point (and then remain 
large afterwards). The expression H5.7|) reduces, in the asymptotic regime, to 

Qf„ w - cot(A^% + const) , (6.4) 

where the const = — 7r/4 can be neglected for the purposes of our argument. Let us denote X^ zs = as + ibs 
with 05, bs real. We then have 

sin as cos as — i sinh 6s cosh 6s , , 

Qm — — TTT ■ (6-5) 

sm as + smh os 

It is clear that if 6s is of order 1 or larger, Qf^ will always be of order 1. Hence 6s will have to be small; we can 
then approximate sinh6s « 6s, cosh 6s ~ 1. For small 6s the expression (|6.5(l will be largest when sin as ~ /3s, 
and the maximum value attained will be ^ l/6s (both the real and imaginary parts will attain values of this 
order). Comparing with the required order of magnitude 1)6. 3|l we find that switching can only occur when 
6s ^ |A^|/A;^. This means that the imaginary part of the wavenumber itself must be tiny: 

zs k^zs 

2 

< fc^as^?^, (6.6) 
Zs 



where we have used that the typical value of the wavenumber is ^ kuA- This condition is not satisfied by our 
model: while (|6.6|l implies that the penetration length Zp ^ (Im {A"^})"^ must be much longer than the SM 
throat length zs, at the end of Section l5.ll we have argued that the penetration length will be much shorter 
than the throat length (and the "self-consistent" penetration length obtained from 1)5.26(1 does indeed come out 
shorter than the SM throat length). Hence we conclude that the mode switching cannot occur. 



7 Enhancing KK mode tunnelling between throats 

In previous sections we have found that, gencrically, energy transfer between the throats will proceed at the rate 
given by the plane wave tunnelling rate of ^14, . In this section we look for ways of enhancing the tunnelling rate 
by modifying the effective potential barrier in the UV region of the RS model. In standard RS the tunnelling 
probability of a mode of mass m (with mL ^ 1) is suppressed as (21 

PRS-{7yiL)\ (7.1) 



Throughout this section we will, for convenience, use curvature radii L instead of curvature scales fc; they are 
straightforwardly related by 

L^\. (7.2) 

Since we are interested in plane wave tunnelling rates, we will work with the plane wave boundary conditions 
that are commonly used in tunnelling calculations, i.e., the SM side wave function will be an asymptotically 
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purely outgoing wave \J mzH^ {mz) . On the A side we will have a mix of incoming and outgoing waves (with 
a nonzero net incoming flux) , V mzH2 {mz) and \/mz7J^ {mz) respectively, where 

H^{x)^J,{x)±iY,{x) (7.3) 

are the usual Hankel functions. Unlike the Z2 boundary conditions employed in previous sections, the tunnelling 
ones do not explicitly depend on the location of the end-branes. We will not take into account the effects of 
decay here; the decay effects can be obtained by a straightforward adaptation of the results of previous sections. 

Our aim of lowering the potential barrier will be achieved by replacing the (single) UV brane with ten- 
sion l|2.5|l by a certain number, N , of branes with smaller tension 

T' = T/N . (7.4) 

Such "generalized Randall-Sundrum" setups were considered, e.g., in ^J^]. If N is small, say 2, tunnelling can 
be speeded up by resonant effects akin to the resonant tunnelling phenomenon known from ordinary quantum 
mechanics (resonant tunnelling was also discussed by (12.. 13 ); we give an explicit example of such a setup in 
Section [7TS1 On the other hand, a large number of appropriately placed branes will lower the effective potential 
barrier and hence enhance tunnelling, as we show in Section f7.5l All cases will have in common the fact that the 
curvatures on the two sides of any particular membrane differ; therefore we will start with a general discussion 
of the matching conditions for such situations. 



7.1 Many-brane Randall-Sundrum 

Let us assume that our system consists of N AdS regions bounded by A^-|- 1 branes located aX z = Zi,i — 0, ..,N . 
The brane at Zq will be taken to be the annihilating brane that are responsible for inflation, while the one at 
^AT+i will be the SM brane. In each region between the branes the metric is the standard AdS metric (|2.2|) . but 
the warp factor will have different coefficients in each region, namely 

a.(z) = -^. (7.5) 
z + Li 

The variables Li are (up to sign) the curvature radii of each region. The sign of Li will be positive (negative) 
if the warp factor decreases (increases) with increasing z. The shifts Li are chosen such that the warp factor 
is continuous at every membrane. The actual values of Li depend on a choice of the origin of the z coordinate 
and the normalization of the warp factor. 

The curvature jump conditions at each of the branes require that the reduced tension Ti = Ti/M^ of brane 
i and the curvature parameters Li and Li+i on both sides of the membrane satisfy (cf. (|2.5(l ^ 

f, = 12 (-^ 1) . (7.6) 



-Ji+l 



L, 



This condition is valid both at the interior branes as well as at the inflating and SM branes, the only difference 
being that at the latter two we also impose Z2 boundary conditions on the background. 

The mode equation in each region has the same form as in the single UV brane case, namely H2.8|l . but with 
the potential V{z) now given by 

Each brane will contribute a term of the form 

V}'^ = f,a{z^)S{z - z,) (7.8) 
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into the effective potential as well. 

Away from the branes in each region the KK wave functions ijjm{z) are given by a general combination of 
the Bessel functions J2{'fTiZi{z)) and Y2{m'zi{z)), where 



Zi{z) = \z + Li\ . 
At each membrane the wave function ipm must be continuous, 

1pm{z "> Z^) = Vm(2 Z + ) , 

and its derivative must obey the jump condition 

The zero mode wave function is a direct generalization of 13. 4|) : in i-th region it is 

^o(^) = ^ 



(2,(z))3/2 ' 



(7.9) 
(7.10) 

(7.11) 
(7.12) 



where the normalization constants Ni are chosen such that the wave function is continuous across every mem- 
brane. ^-'^ In fact, since the warp factor a has the same functional form and coefficients that do make it continuous, 
we can also write 

M^) = NoCT^'/^iz), (7.13) 
where the normalization constant A^o is now the same for all regions. 



7.2 Throats with different curvature radii 

An interesting generalization of the models of [HIE] is to allow the curvature radii of the two throats to differ. 
The radii can in principle arbitrary as long as the tension of the brane that separates them obeys the jump 
condition (|7.6|l . Let us first look at the plane wave tunnelling rate through a one-brane potential barrier that 
separates two regions with curvature radii and , respectively. Normalizing the warp factor to 1 at the 



UV brane (z = 0) leads to 
are, respectively, 



L^, — (see H7.5|l 'l. The KK mode wavefunctions on the A and SM sides 



■\/ mzA [-AH2 (mzA) + BH2 (itlza)] 
^/mzsCH^imzs) , 



(7.14) 
(7.15) 



where 



ZA,s=\z + L^'^\. (7.16) 

As usual, on the A side we have a superposition of an incoming and a reflected wave, whereas on the SM side 
only the transmitted component is present. Continuity and first derivative jump at the brane, H7.10|l and (|7.6|l . 
lead to 

C _ flAH2{mL^)H+{mL^) - H^; {mL'^)H+ {mL^} 
A ~ ' 



(7.17) 



H+{mLS)H+{mL^) + H+{mL^)H+{mL'S) ' 

Let us first assume that the curvature radii arc such that mL^ <^ 1, mL^ <C 1. Then, using the small- argument 
expansions of Bessel functions (|A.ip , we find 



to2(L^)3/2(^S)1/2^ 



(7.18) 



'^The proper jump in the first derivative then follows automatically. 
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Next, let us assume that the SM curvature radius is large, such that mL ^ 1. The Hankel functions can 
then be characterized as H^{mL^) ^ 6^"""^^/%/ mL^ , and we find 



In particular, this formula also applies when oo, that is when the SM side of the UV membrane is flat, 

and so characterizes the tunnelling of KK modes out of a single AdSs throat. Compared to the equal radii case 
= , the tunnelling rate P — \C/A\^ is clearly enhanced by a factor of l/(mL"*). 

Having significantly different radii in the two throats has other interesting consequences: for example, it can 
relax the constraint H3.29|l on the relative sizes of the warp factors at the A and SM branes - the derivation 
of the constraint from H3.28|l (which remains valid even in the case of different curvature radii) assumed equal 
curvature radii on both sides. Conversely, if the curvature radius of the SM throat were much larger than 
L^, the SM warp factor 175 = /{L^ + zs) could be much larger than the bound (|3.29(l while allowing KK 
modes to tunnel from the A throat to the SM one. Phenomenologically, however, such a setup may not be 
desirable because of the presence of a large hierarchy between the curvature scales in two throats; it may also 
be difficult to realize such a setup as a string theory compactification. 

7.3 Resonant tunnelling through a gravity box 

One possibility to speed up the energy transfer between throats is to use an analog of the resonant tunnelling 
effect from standard quantum mechanics: split the single potential barrier (the UV region) in the RS geometry 
into two, with a flat region (a well with zero potential in the effective Schrodinger equation) inbetween . If 
the frequency of the tunnelling wave is tuned to the size of the of the well in the middle, the tunnelling rate can 
approach one. As argued in |12l I13| who analyzed a similar setup, in a realistic string theory compactiflcation 
one would expect that the two throats are joined by a bulk region that is modelled by the flat interval. 

Hence we consider a setup of the type described in |15| . but with the tunnelling boundary conditions. The 
geometry consists of two AdS throats, whose curvatures we will take to be the same for simplicity, separated 
by a flat region. The curvature radius of the throats will be denoted by L, while the length of the box in the 
middle will be 21, with the branes sitting at z — ±1. The effective potential is plotted in Fign(a). The warp 
factor is constant in the central region and will be normalized to 1 there (so the coordinate and proper length 
of the middle box coincide). On the A side (region 1) the wave function of a KK mode with mass m is 



C 
A 



(7.19) 



ijl^\z) = ^/mzi{z) [AH+ {mzi{z)) + BH^ {mzi{z 



))], 



(7.20) 




C coa{mz) + Dsm{mz) 



(7.21) 



and on the SM side (region 3) we have 




y/mz3{z)EH^ [mzsiz)) . 



(7.22) 



In the notation of Section 17.11 the parameters Li are 



L3 



Li = L, L2 = L2 



00 , 



(7.23) 



while the shifts are (from requiring a{z — ±1) = 1) 



L3 



L + l. 



(7.24) 
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It is straightforward, if somewhat tedious, to solve the matching conditions H7.10|l and H7.11|l . To shorten the 
resuh, let us denote 

i EE tan mZ, (7.25) 

and omit the arguments of all Hankel functions, since they are all evaluated at rriL. The tunnelling amplitude 
E/B is then 



E_l {t+\){H^Ht -H^H, 



(7.26) 



As a check, when i ^ 0, so 1/t dominates over everything else, the above tunnelling amplitude reduces to the 
bulk- less expression given in |14j . 

The tunnelling resonance occurs when the denominator becomes small, i.e., when 

If we concentrate on KK modes whose mass m is small compared to the inverse curvature radius, mL <ti 1, we 
can use the asymptotic forms (|A.ip of the Bessel functions, and we find that the denominator becomes small 
when 

i - t - mi < 1 , (7.28) 

that is when i w 1. This is not surprising ~ it means that to get an amplification of the tunnelling rate, the 
size of the central box must be roughly the same as the wavelength of the tunnelling particle. As a numerical 
illustration, we plot the potential and the tunnelling rate of this setup in Fig. ^ The parameters used in the 
plot are / = 5, L = 1. 
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Figure 1: Two branes with a flat region between them. Graph (a) plots the potential V{z), while (b) shows 
the tunnelling rate P as a function of the KK mode mass m (solid line) together with the reference tunnelling 
rate of the standard double RS geometry (dotted line). The brane separation is 21 = 10. Tunnelling is greatly 
enhanced for masses that fall within the resonant tunnelling peaks. Note also the suppression of the tunnelling 
rate for very light modes. 

This is bad news for phenomenology. According to the discussion around Eq. (|2.21() we would expect the 
masses of the produced KK modes to be smaller than the annihilation scale Ma, which in turn is expected to 
be a few orders of magnitude lower than the inverse curvature radius L of the throats. Hence, the flat region 
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would have to be unnaturally long as compared to the characteristic length of the throats, and presumably 
as compared to the radius of the internal manifold in a string theory compactification; in other words, the 
hierarchy problem that the warping was supposed solve would reappear in a different guise. 
Lastly, one can consider the regime of t small but not too small; more precisely, 

l>t> ' ' r^mL. 7.29 
Y2[mL) 

In this regime the tunnelling amplitude behaves roughly as 

^-—imL)\ (7.30) 

meaning that the tunnelling rate with bulk present is suppressed by an additional factor of (mL/t)^ compared 
to the bulk- less tunnelling rate. The extra suppression is also evident in the plot^b). 



7.4 Comments on the WKB approximation 

Incidentally, our results H7.19|l and (|7.30|l provide a cautionary note on the WKB methods used, e.g., by jl2lli;-{j . 
In the WKB approximation, when the tunnelling amplitudes are small, the tunnelling amplitude Qtt between 
two throats is roughly given by the square of the rate @tb at which particles can tunnel from one throat into 
(flat) bulk, 

Qtt-e^. (7.31) 
It is clear that such formulae must be used with caution: for example, the tunnelling amplitude between two 
RS throats, calculated in 14 , is Qtt ~ [mLY , while the throat-to-bulk tunnelling amplitude H7.19|l is Qtb ~ 
(toL)'^/^. Further, when the two throats are separated by a flat bulk region, the tunnelling amplitude l|7.29(l does 
formally exhibit the {mL)^ behaviour that one would expect from Qfb, but it also depends on the length of the 
bulk region represented by the parameter t (and the dependence is very different from the WKB approximation 
as used in [12 El). 



7.5 Lowering the potential barrier by multiple branes 

Another way of speeding up the tunnelling process is to lower the potential barrier in the UV region by replacing 
the single UV brane by a large number of branes of lower tension that are also spatially separated along the 
z direction. The basic idea is to utilize the fact that the (squared) warp factor (jf{z) and the potential V{z), 
given respectively in (|7.5|l and H7.7|l . differ by a factor of the (squared) curvature radius L^. By arranging the 
branes such that the the warp factor becomes of order 1 only when the curvature radius is also large, we can 
achieve a significant lowering of the effective potential barrier between the throats. 

We have not attempted an exhaustive search for the configuration most favourable for tunnelling. It is clear 
nevertheless that the the modification must happen in the region where the warp factor (and the Randall- 
Sundrum potential) are large, namely around z = 0. In this section we present a toy example of such a 
modification. 

Suppose we arrange N — 2n branes with tension Ti — T/N such that the middle region is (approximately) 
flat. Then the curvature radius of a region i satisfies 

^ = 12fi{i~n). (7.32) 
Li 

As discussed above, l/Lf is precisely the factor (up to 15/8) by which the potential and the warp factor differ: 

^•w = fir 
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The last expression makes the suppression of the potential more apparent: the warp factor is largest in the middle 
(flat) region, but the potential there is suppressed by the {i — nY /N'^ factor. Because the shifts Li entering 
the warp factor are determined recursively, it is not possible to derive a simple formula for the maximum of the 
potential (|7.33|) . so we give a few simple examples for illustration. 

In our example, we fix the warp factor value at the maximum to be 1. This choice ensures that the warping 
at the "top" of the throats is the same for all setups — while arbitrary, this has the physical consequence of 
measuring masses at the top of the throats in the same units. Further, since we consider a many-brane scenario 
to be a near-the-top modification of the "standard" model of ^3 , we fix the curvature and warp factors at 
the bottom of the throats to be the same in all cases. 

In Fig. 12 we show potential as a function of the coordinate z and the tunnelling rate as a function of the mass 
m of the tunnelling KK mode for the benchmark doubled RS (i.e., constant-curvature, single central brane) 
setup of jni[I3|. As expected, the tunnelling rate becomes substantial when the KK mode mass m? becomes 
the same as the potential at the top of the barrier. 
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Figure 2: The single central brane setup with constant curvature. Graph (a) plots the potential V{z)^ while (b) 
shows the tunnelling rate as a function of the KK mode mass m. 

Now let us look at a setup where the central brane is replaced by 100 equidistant branes (50 on each side 
of z = 0) with tension equal 1/100 of the tension of the single brane. The potential and the tunnelling rate are 
shown in Fig.|31 In Fig. ^ we zoom in on the central region to show the potential in more detail. The potential 
barrier is significantly lower and the tunnelling rate for low-lying KK modes is correspondingly enhanced. 

Are the implications of this toy model applicable to an honest string theory compactification? Recent 
work |12[I13| studying string theory throats seem to answer this question in the negative. We would like to take 
a more cautious position: it seems clear that the tunnelling rates depend to a large degree on the details of the 
gluing of the throat to the bulk; these details that are not yet under sufficient control. In fact, the effective 
potentials for the KK modes in ^1^] and our potential plotted in Fig. O are qualitatively similar, but their 
quantitative details are apparently sufficiently different to cause a major difference in tunnelling rates. 

8 Conclusions 

We have presented an exhaustive discussion of the properties of the KK modes that are responsible for the 
energy transfer from the A throat to the SM throat. We have found that tunnelling in the conventional sense 
will proceed only if the SM throat is substantially longer than the A throat, namely only for zg ^ ZAikzA^. 
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Figure 3: The multiple brane setup. Graph (a) plots the potential V{z), while (b) shows the tunnelling rate as 
a function of the KK mode mass m. 

Even when the A throat KK modes do not tunnel in the conventional sense, their wavefunction has a tail in 
the SM throat that is big enough to induce a large decay rate into 2 SM throat KK modes. Naively this rate 
should increase with the length of the SM throat and, for a long SM throat, will be large enough to make most 
of the annihilation KK mode energy go into SM throat KK modes. However, a more careful analysis taking 
into account the effects of wavefunction decay shows that the decay rate into two SM throat KK modes cannot 
get larger than a certain maximum. The effect can be understood intuitively from the fact that a large decay 
rate would shorten the penetration length of the A throat KK mode into the SM throat, which in turns leads 
to a suppression of the overlap of its wavefunction with the SM throat KK modes, and that in turn suppresses 
the decay rate. We have shown that the physically realized decay rate is effectively given by the plane wave 
tunnelling rate. Nevertheless, the naive decay rate is not completely meaningless: we show that it is turn small 
enough that it does not induce a further slowing down of the tunnelling/decay rates via the complex frequency 
effect discussed in (10) . 

An additional mechanism of energy transfer we have investigated is based on level repulsion and can be 
pictured as throat switching by KK modes when the spectrum of the SM throat is gradually lowered by the 
relaxation of the SM throat. We found that the presence of complex parts in the mode frequencies destroys the 
level repulsion effect, because levels can now avoid each other in the complex plane; hence, the energy transfer 
by switching cannot occur. 

We have then proposed a simple modification of the UV region of the geometry that, at least in the phe- 
nomenological 5-dimensional model, leads to a drastic enhancement of tunnelling rates between the two throats. 
The modification consists of replacing a single UV brane with an array if branes with lower tension such that 
the curvatures deep in the throats remain the same. This has the effect of lowering the potential in the effective 
Schrodinger equation that the KK modes have to tunnel through, leading to an enhancement of their tunnelling 
rates. It is not clear, especially in the light of the recent studies of stringy throats |12[ I13| . whether a similar 
mechanism could be realized in an honest string compactification; we argued that a better knowledge of the 
geometry where the throats are glued to the bulk is necessary to answer this question definitely. 
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Figure 4: Detail of (a) the potential V{z) and (b) the tunnelling rate P{m) shown in Fig. |31 
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A Modes of a two-throat system 

In this appendix we study in detail the modes of the two-throat system described in Section 13.11 
A.l Small mass modes: mzA ^ 1 

Taking the mass to be small allows us to use standard formulae for Bessel functions of small arguments, namely 



We then have 



2 1 

Ji{x) K ]-x 



4 1 

Y2{X) « — - 
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and find 

One should not be surprised that this ratio is not dimensionless; after all, we are comparing wavefunctions 
that have different (namely continuum vs. discrete) normalizations. Physically one cannot talk about a single 
mass value (or mode) out of a continuum; one must always integrate over a certain range of masses to get a 
meaningful result. One can, e.g., integrate the probability density of finding a KK mode of mass m! < m &t 
the A brane and compare it to the probability density of finding the zero mode there, i.e., calculate 

= r^"^;(e(^A))^ 



Using HA.7|) we find 

1 rn? 
4 fc2 

As derived, this expression is only valid for m <ti ^/za. We would like to extend this calculation to modes whose 
mass is of order 1/za, because we expect that the effective brane scale Ma will be roughly of that magnitude. 



A. 2 Medium mass modes: itlza ~ 1 

There is no good approximation formula for the Bessel functions when their arguments are of order 1; the best 
one can do is to say that generically, all of the relevant Bessel functions are also of order 1. Thus one cannot 
say much more about than that it will also generically be of order 1. That is enough, however, to fix by 
the jump condition at the Planck brane (recall that we still assume mk ^ 1, so all Bessel functions evaluated 
at the Planck brane can be replaced by their small-argument approximations). We find that Q^, and therefore 
also the normalization constants and are given by the same expressions as in the small mass case, 
namely (jA.3|l — (|A.5|) . The wavefunction of the massive mode at the A brane is then^^ 

5/2 '^-'-/'^ 

V'™(^a) ^ (A.IO) 

Interestingly enough, the ratio %Ij^{za)/^o{za) is then (using mzA ^ 1) roughly the same as in the small mass 
case, namely 



VoyzA) k 

We show below that this result holds only for truly "generic" values of to; we will see below that at special 
points, the ratio i?^ can grow to be much larger than the above value. 

The result (|A.lip implies that one can integrate up to m ^ za in (|A.8p with result 



Pi^)-[j) ' (A-12) 



that is the non-zero mode amplitude at the A brane is suppressed by roughly the square of the ratio of the 
inflation and string scales (modulo the radius of the 5th dimension). 

Since in this case we do not have factors of order one under control, it is useful to plot the iP:^{za)/iPo{za) 
numerically. The plot can be found in Fig.Efa). 

To see that the validity of our approximations is restricted to "generic points" , let us look at plots of the 
ratio for values of m that are several times I/za- The plot in Fig. |5{b) shows that for certain masses, 

^^We use the ~ sign to denote our ignorance of factors of order 1 (we reserve f» for approximations that hold up to small 
corrections). 
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Figure 5: Plot (a) shows the ratio i?^ = ■iP„^{za)/iPo{za), as a function of m (in units of 1/za)- For small 
m the ratio behaves as ~ y/m/k. In this plot k = 1, za = 10^. In (b) we plot the logarithm of the ratio 
— ip^{zA)/i'o{zA), as a function of m (in units of 1/za), for a larger range of m. This plot is logarithmic, 
so the KK mode wavefunctions are suppressed with respect to the zero mode whenever the curve is below the 
horizontal axis. Note the presence of spikes that actually turn out to dominate the integrated ratio p{m). 



the wave function of the corresponding modes on the A brane is parametrically amplified. One might suspect 
that these spikes in i?^ appear when either Q^^ diverges (at the zeros of Ji(toz^)), or the denominator in the 
expression (|2.18|l for becomes zero, and R^, would diverge. In fact, it turns out that neither of these special 
points can cause the spike; rather, the spikes appear at values for which approaches zero, and R^, remains 
bounded (but large) . For m <^k one can use the approximate values ljA.l|) of the Bessel functions to show that 
at the points where Qf^ = 0, the KK mode wave functions attain values ~ k^ jvr? , and the ratio R^ grows to 
R^ - {k/my/^/^. 

Physically, the interesting question is for what values of m does the integrated probability density ratio p(rn) 
of (|A.8p reach values ^ 1. The presence of the spikes makes a numerical integration difficult, so we provide 
estimates that will turn out to be good enough. We start with generic values of m and consider the contribution 
of the spikes separately below. 

For a values of the KK mode mass m away from the spikes, the approximate form (jA.10|l holds as long as 
m < k. Keeping za in the expressions explicitly, we arrive at 

R^-^y^imzA)\ (A. 13) 
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and 



p{m) 



m 
fc2' 



{itiza) 



(A.14) 



This result is counterintuitive: it says that for the KK modes to dominate over the zero mode, their mass must 
be so large that they are already far into their asymptotic regime, m 3> 1/za- However, in the asymptotic 
regime we would expect the KK modes to dominate over the zero mode, because the KK modes approach plane 
waves with constant magnitude, whereas the zero mode falls off as a power of z. We seem to be missing some 
KK mode contributions; this discrepancy is a hint that the spikes contribute significantly, and actually dominate 
p{m). Let us now show that. 

As we have already mentioned, the spikes appear at masses m such that Q^^ — 0. We remind the reader 
that we still assume fh k. From H2.16|l we find that the corresponding must be 



Yi{m/k) 
Ji{rfi/k) 



^ > 1 



(A.15) 



This equation together with (|2.14|l determines the values of m where the spikes occur; the first few spikes will 
be at TO 1/ za <SC fc in accord with our assumption, ft is also worth noting that the spike occurs when Ji{mzA) 



is small but not zero, Ji{mzA) ^ ^/Qr, 



From H2.f 7|l we find 



(A.16) 



and from (|2.18|l . taking into account (|A.15(I 

Ni^Ng. (A. 17) 

To estimate the contribution to p{m) from a spike, we will find its width Sm at half maximum. The relevant 
quantity here is R'^{m) = iI}^{za)/'4>'q{za), so half maximum is the point where the KK mode wavefunction 
squared "0^ drops to half its peak size. The drop can come from two sources (see (I2.10f) '): either or can 
halve. Let us for now assume that varies faster, find the width 5m and then show that the variation in 
due to 5m is smaller. 

Hence, we are looking for 5m such that for to = m + 5m, (iV^)^ = (iV^)^/2. This implies (iV;fJ^ = (N^)'^ /2 
and hence 

Ql^l. (A. 18) 

Keeping only the leading terms in (|2.16(l (recall again m/k ^ 1) we find that must change from the 
-2Yi{m/k)/Ji{m/k) of (TOsjl to 



2Yi(TO/fc) - Ji{m/k) 
Ji {m/k) 



1 



(A.19) 



that is 5Qm ~ — 1. Recall that is large because Ji(toza) is near its zero; thus 5Q^ will come from the 
change in JiimzA) caused by 5m, while Yi{mzA) remains approximately constant. Denoting the change in 
Ji{mzA) by 5Ji, we require 



leading to 



Taking into account Yi{mzA) 



YijmzA) _ / 
Ji{mzA) + 5Ji 

5Jl K 

1 we have Ji (fhzA ) 



1 



Yi{mzA) 
Ji{mzA) 



1 



Jl{mzA) 
Yi(toza) 



(A.20) 



(A.21) 



leading to 



(A.22) 
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The last step is to write SJi^ SmzAJ^ijnzA) ^ SmzA (the derivative of Ji near zero is ~ 1) giving 



Thus each spike will contribute roughly 



IT) ^ 1 



5p ~ 6mR^{m) 

^ (A.24) 



kzA 

into the integrated wave function ratio p{m). If we concentrate on the first few spikes where rnzA ^ 1, we can 
also express 5p as 



op ^ — 

m 



2 

2 



{kzA 

^A 

The last form is particularly interesting - this is precisely the enhancement factor used in the literature [5ll51llO| 



a^^ . (A.25) 



B (Non-) Tunnelling between two wells in quantum mechanics 



In Section 13.31 we have stumbled on a surprising fact: particles will not timnel from the shorter into the longer 
throat unless the longer throat is substantially longer; in the AdS throat case we found that ratio of the throat 
coordinate lengths, zs/za, had to obey 



ZA 



> 



k 



(B.l) 



A natural question of interpretation arises: in general, given a potential barrier between two wells, how much 
longer must the target well be for particles to tunnel from the source well? The answer is quite important 
in the context of warped reheating where the potential in the UV region presumably differs from the simple 
doubled RS model; if such a modification changes the above inequality (or, equivalently, 13.28|l 'l. the conclusions 
of Section rOl would be changed as well. 

To gain some insight into the problem, let us consider a simple 1-dimensional quantum mechanical problem 
of two wells separated by a 5- function potential barrier located at z = 0, 



Uiz) = UoSiz) 



(B.2) 



In analogy with the two throat model we call one of the wells, for z < 0, the "annihilation" well and make it of 
length ZA, the well on the positive side will be called the Standard Model (SM) well and will have length zs- 
Note that za and zs correspond to their respective namesakes in the two-throat model which are coordinate 
lengths, not the physical lengths y '^Inz. Hence we will be studying the Schrodinger equation 



2dz 



n— + UoS{z) 



ij{z) = ^E^^,{z) 



with the boundary conditions 
At z = we require continuity, 



iP(z = -za) = i^{z = Zs) ^0. 

^p{z = 0+) = v(z = o-). 



(B.3) 

(B.4) 
(B.5) 



32 



and the appropriate jump in the first derivative, 

2 W) 



Uo . (B.6) 



Analogously to the AdS case, let us denote the wave function corresponding to energy E hy tpE] for clarity, we 
will write ipEiz) = ipsi^) z < and iPe{z) = 4'e{z) for z > 0. The Schrodinger equation (|B.3I) away from 
z = has the solution (taking into account the boundary conditions (|B.4|) ') 

V'bW = Ni sm[Eiz + za)] 

^l(z) = Nism[E{z-zs)]. (B.7) 
The matching conditions at z = then imply 

sin EzA = - N§ sin Ezs, (B.8) 

cot EzA + cot Ezs = -2^- (B.9) 

E 

The last equation, (|B.9|) , is the one that determines the spectrum. The "barrier" regime of parameters is when 
Uo ^ E; on the other hand, if Uq <C E, the mode wavefunctions will not feel the presence of the i5 function 
potential appreciably. Let us therefore concentrate on the barrier regime Uq ^ E. Then the eigenvalues of E 
will be those for which either cot Eza is large and negative (we call these the A side modes) or where cot Ezs 
is large and negative (these modes will be called the SM side modes). We do not consider the case where the 
well lengths are tuned and both cot terms become large at the same time. 

Let us first look at the A side modes. We have cotEzA ~ —2Uq/E ^ 1 so sinEzA ~ ~E/{2Uo) while 
cos Ezs ^ 1. Then the continuity condition (|B.5|) implies 

^E-^E^^- (B.IO) 

We then find that the probabilities and Pj of finding the particle on the I and SM side, respectively, obey 

p| _ {Nirzs 



Pi {NirzA 

zs f E 



(B.ll) 



ZA \2Uo 

We remind the reader that the corresponding expression in the two throat RS setup was 

f 5 (f )* ^ 

The two formulae seem to be very different; in particular, the powers of E and m (which are analogs of each 
other) differ. 

This discrepancy brings us to a question we should have answered right in the beginning: how do we 
compare the two setups, or, more precisely, what parameter values in the square well model should give results 
comparable to a given set of parameters in the two throat RS setup? A naive guess might be that the integral of 
the potential in both cases should be comparable; however, there is no good argument saying that the integral 
of the potential is a physically meaningful quantity. Hence we propose to set the model parameters such that 



^3 In fact, if one considers a slowly varying potential W, one can approximate the wavefunction as exp[\/VF — E'^z], so the relevant 
quantity describing the suppression of a wave function by the barrier would appear to be exp[J dz\/W — E'^] , i.e., the integral of 
the square root of the potential. Such an integral does not make sense for the 5 function potential wc arc considering; of course, 
the approximation of a slowly varying potential is obviously invalid as well. 



33 



the tunnelling rate of plane waves (which is a mathematically well-defined and physically meaningful quantity) 
are the same for the two barriers. 

In fact, the tunnelling rate Prs of the doubled RS model is known and is given precisely by the factor (m/k)^ 
entering (|B.12p . Let us therefore calculate the tunnelling rate for the (5- function barrier. We are considering 
plane waves and hence effectively infinite za and Z5; the problem now has a continuous spectrum. For the 
tunnelling calculation we simply require the wavefunction on the incoming (A) side to be a combination of the 
incoming and reflected plane wave, 



while on the outgoing (SM) side the wavefunction should have just the transmitted component. 
Continuity at 2: = implies 

A + B^C 

while the first derivative jump condition leads to 

Extracting the tunnelling rate P,^ = C"^ /A^ is simple and the result is 

^0 



B.13) 

B.14) 
B.15) 

B.16) 
B.17) 



Hence the tunnelling rate Pgq is, up to a factor of 4, precisely the factor entering ljB.ll|) (we disregarded such 
factors in the derivation of IjB.llj) anyway). 



C Decay rates of KK modes 

We present an effective 4-dimensional calculation of the decay rates of KK modes into two particles that can be 
either gravitons or lower lying KK modes (or one of each). The effective 4-dimensional 3-point couplings C for 
these processes are calculated in the main text. The calculation we do here is in principle a standard textbook 
one with the added complication of having a tower of KK modes accessible as decay products. 
The effective 4-dimensional interaction is 

Sint^ J d^xC(t>Q{x)(j)i{x)(j)2{x) (C.l) 

where we take (j)Q to be the decaying mode with mass toq a-nd 0i, 02 are the products with masses rni,m2 
(which need not be nonzero). Note that the coupling "constant" C can depend on the momenta of the particles. 
Since the decaying mode is massive, we can go into its rest frame. Let us denote the particles' energies by 
ojq = mo, (^i,uj2 and their space momenta by po — 0,pi,p2- The decay amplitude is simply ^ and the decay rate 
^mi,m2 is (we omit all factors of 27r) 

Fmi.m^ / d^Pl(fip2S^iPl +P2)S{mo - UJl- OJ2) . (C.2) 

J mQUJiuj2 

When comparing the decay rates into KK modes and into gravitons, one must perform a sum over all accessible 
KK modes to obtain the total decay rate Vkk into KK modes. 



^2KK = 2^ mi ,7712 



mo pm—mi 

dnii / dm2Tmi,7n2 , (C.3) 
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where in the second Hne we have approximated the sum over a dense discrete spectrum by an integral. If one 
of the final states is a graviton, there is only one integral (say over rrii): 

/•mo 

^KK,g ~ Zs j dmiT,n^.m2=0 (C.4) 

Jo 

Evaluating Tmi,m2 is simple. The integral over p2 is trivial; the energy S function then forces the magnitude 
of pi to be 

^ _ (mg - m,l - ml) - Amfrnj 
Taking account of the S function then gives 

where if it depends on the momenta, is evaluated with pi = —p2 given by (|C.5p (by rotational invariance ( 
cannot depend on the direction of pi). The two main cases of interest in Section 0] are when C is independent 
of momenta and when it contains a factor of po ■ Pi = mo-\/ mf + p\. (One could also contemplate pi ■ p2', this 
coupling should be comparable to po ■ Pi-) Let us treat both cases together by writing ^ as 

C^CiPo-Pif (C.7) 

with = 0, 1. We can then write 

^2KK ^ zg I dmi I dm2C,'^ml^~'^{m\+p{Y 
Jo Jo 



where we have denoted 



z^mo'^C^ / dvi dv2Wi{vl + wi{vi,V2f)'^ , (C.8) 
Jo Jo 



i.e., the momentum pf in units of mg. The phase space integral is now written as a dimensionless integral that 
will simply give a number that we expect to not differ significantly from 1; we duly disregard it and obtain 

r2KK zlmt^^'C' ■ (C.IO) 

The total decay rate into one graviton and one KK mode can be treated similarly; the result is simply one less 
factor of zsmo because of the missing integration over m2: 

TKK,g- zsm^f'C^ (C.ll) 



D Inclusion of a sink in a two-well system 

In this Appendix we discuss two- well systems in which the wavefunction decays with rates that are different on 
each side of the barrier. We will first gather basic facts for a single well. 
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D.l A potential well with a sink 



Let us consider quantum mechanics of a square infinite well of length zq that has a constant imaginary term 
in the potential. We write the wavefunction as '^{z,t) — T{t)ijj{z) with T{t) — exp{iujt) and allow ui to be 
complex. The Schrodinger equation for tp is then 



2 



■Sq 



t/j = —Uj'^lp J 



and the boundary conditions are the standard 

4>{z = 0) = = i) = 0. 

The solution is 



ip{z) = A^sinAz 



where 



A — \J llP- — 2iSo — — ■ 
ij 



(D.l) 

(D.2) 

(D.3) 
(D.4) 



We hence find that unlike in the usual case, the wavenumber A and the frequency differ and while the wavenumber 
is real (this is dictated by the boundary conditions), the frequency is complex. 



(D.5) 



and the wavefunction can either decay or grow exponentially depending on which branch of the square root 
one takes. Of course, if the iS^ term is supposed to represent a sink, one must choose the sign for which the 
wavefunction decays. 



D.2 (5- function barrier with a sink 

Let us now consider the toy model of a i5-function potential barrier, where we add a sink term that has a different 
magnitude on each side of the barrier, S = = const for z < and S — — const for z > Q. We again 
write the full time-dependent wavefunction as ^(z,i) — T(t)'ijj{z) with T{t) ~ exp(zct;i), where lu is in general 
complex and is the same on both sides. The boundary conditions ■ip{z — —za) = ipiz = zs) = together with 
the Schrodinger equation IjD.ip (with the appropriate sink on each side) imply 

^{z) = 7V^sinA^(z + ZA), z<0 

^{z) = 7V^sinA^(z-Z5), z>0 (D.6) 
with the wavenumbers given in terms of lu as 

A^-^ = Vt^2 _ 2iS^,s . (D.7) 

The matching conditions at z = are exactly the same as in the sink-less case discussed in Appendix ^ and 
can be written as 

iV-^sinA-^ZA = -TV'^sinA^zs, (D.8) 
A^cotA^ZA + A^cotA^zs = -2Ua . (D.9) 

Let us again assume we are in the "barrier" regime Uq 3> 1/za > ^/zs- To satisfy (jP.Qp . one of the cot functions 
must be large (we assume the wells are not tuned, so at most one of the cot's can be large for any value of lu). 
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As a function of a complex variable x, cot x is large only near x„ = nn. Let us assume that it is the A side cot 
that is large, so we find that the A side modes will have wave numbers 



X-^^. (D.IO) 



In this case the corresponding A'^ will be complex, 



X'n^\ -^~2t{SS~SA). (D.ll) 



The presence of a complex term in the jump condition (|D.9|) implies that will also have a small imaginary 
part and hence that even if — 0, the wavefunction of an A mode will slowly decay via "seepage" and 
subsequent decay on the SM side. The imaginary part of A^ can be estimated as 

zsiS^-S'^) for|5^-5^|< 



Ur, Z , 



I v/^^ - Sa for > rv 



SO even for large it is suppressed by 1/{UoZa)^ , which for low-lying modes (n ~ 1) roughly equals the 
tunnelling probability P^q given in IIB.17|) . 

Let us now look at the case when the sink is only turned on the SM side, i.e., — 0, while S^' ^ and 
comparable to ^/z^. While A^ is (nearly) real, Af = a + z6 has a sizable imaginary part; this means that the 
wavefunction on the SM side, 

V'(z) =sinA^(z-Z5) (D.13) 

will actually be largest near z = and will decrease in magnitude roughly as exp6(z — zs). This is what one 
would expect intuitively: the wavefunction has roughly constant magnitude on the A side; on the SM side very 
close to the barrier the magnitude will be determined by continuity (jP.Sjl and so will be not affected by the 
sink very much, but further away from the barrier the wavefunction decays roughly exponentially (actually as 
a sum of a cosh and sinh term) due to the presence of the sink. As a consequence, the wavefunction remains 
normalizable to 1 even in the limit zs c>o (while keeping za constant), rather then becoming plane wave 
normalizable as in the case without a sink: continuity HU.8|I implies 

|7V^|2 = liV^psinA-^z^i^ 1 



cosh2&Z5 ' 

while the normalization integral on the SM side is 



I sin(a + i5)z5p 

iTV^HsinA^ZAp , (D.14) 



/•zs 

tS 



^ dz| sin(a + i6)(z - zs)| 

"'0 
/■zs 

= dz [cosh 2b{z — zs) + cos 2a(z — zs)] 

Jo 



« ^sinh26zs (D.15) 

so the product \N^\^I^ remains finite for zg oo. Effectively the particle penetrates the SM side only to 
distance z ~ 1/6 (that penetration is of course suppressed by the usual tunnelling probability - the factor 
sinA^z^ is small). The magnitude of b depends of course on we remind the reader that a + ib — X ~ 
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E Sinks and reduced decay rates 



In this Appendix we derive the values of the sinks S^, in a two-throat RS system that reproduce given decay 
rates F^^jF'^. The sinks are an effective description of the decay; in particular, they are different for different 
frequencies. Indeed, if the sinks were frequency-independent, the decay rates for heavy modes with mass m ^ S 
would be given by F ~ S/m (where the frequency lo = m + iV) and hence would decrease with increasing m, 
contrary to what one expects, and obtains, for the decay of KK modes in our model. Apriori it is not completely 
clear that the sinks have to be the same for two modes of the same frequency that are localized in different 
throats; however, if the sinks were different, one runs into potential inconsistencies with some modes simply 
disappearing from the spectrum when one of the throats changes length, so it appears that the sinks cannot 
depend on whether the particular mode is localized on one side or the other. In fact, this conclusion can be 
verified independently by an explicit calculation based on the continuity equation of the probability density; we 
will not describe the calculation in detail here. 

Let us take a particular mass and denote 6S = ~ S^; then we have (A'^)^ = (A"^)^ -l- 2iSS. Given the 
decay rates F"^ and F'^ of a particular mode in the A and SM throats, respectively, we will find useful to define 
reduced decay rates Fq ,Fq that only depend on the properties of the particles the mode decays into, while the 
dependence on the size of the wavefunction of the decaying mode as well as on the length of the throat is scaled 
out: 

F^ = r^|(iv^f(i + (g^)2)|%^, 
F^ = ro^|(iv^)2(i + (g^)^)|^5. 

For example, for an A throat mode \ iN^)^[l + (Q^)2]|za 1, so F^ - F^, but |(A^^)2[1 -t- (Q^)^]ks ^ cri/c^^ 
and hence F"^ ~ T^lX/kl'^zs/ za- In other words, the reduced decay rates encode properties of the throat and 
modes into which the decaying mode decays into. Regarding the two-throat system as two weakly coupled 
single-throat systems leads us to expect that the reduced decay rates should determine, up to small corrections, 
the corresponding sinks (and the corrections should vanish in the limit of an infinite barrier) . 

Let us start by considering an A throat KK mode. The aim will be to find the imaginary part of its 
wavenumbers induced by the presence of a barrier and a particular set of sinks; this will allow us to express the 
(known) decay rate of the mode in terms of the two (unknown sinks) . Repeating the same procedure for an SM 
throat mode of the same mass (that has a different decay rate but the same sinks) gives us the second equation 
for the sinks; this system can be then solved to find the sinks. We will assume that 6S is at least somewhat 
smaller than m^, and we will work to lowest non-trivial order in perturbation theory in m/k. 

For an A throat KK mode, the matching condition (|5.8I) implies 

\SvSvS I \AvAvS 
Q^n XAjAyS ' ^^-^^ 

where we have denoted = Yi{X^ /k) etc. Using the small-argument approximations (|A.1|I we have 

^™*"(W(^ + W)' ^^'^^ 
so is large. Using the general formula (|2.14() and the asymptotic form of the Bessel functions we have 

Qi « cot{X^iA ~ 3^/4) ; (E.3) 

\Qm\ ^ 1 tl^™ implies that 



sin(A^iA - 37r/4) 



{X^f 1 

^ + (PT)2 



(A^l / 



i6S 



P 2 V (A^) 
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This in turn implies that 



37r/4: 



(A^l 
B 2 



1 



i5S 



(E.4) 



If the two throats were decoupled, we would have X^^^za — 37r/4 « nvr; let us denote the imaginary part of A"^ 
induced by the presence of the SM throat with its sink by 5\"^ = Im {A'^ — A^^^}. Eq. l|E.4p can then be written 
as 



6X' 



SS 



(E.5) 



The derivation here used the assumption cos(A'^z^ — 37r/4) w 1; this is true whenever the imaginary part of 
X^Za is small. Under our assumptions we indeed have 5X^za ^ 1, so the calculation is consistent. 

The imaginary part of the wavenumber determines the decay rate Ta of the mode via (recall that (m + 



iVAl^y 



(A 



A\2 



Ta ~ 25X^ 



2S^ 



(E.6) 



In general, the decay rate of the (A throat) mode Ta is the sum of the decay rate in the A throat and 
the decay rate in the the SM throat, Ta = F^ + F^. We would like to express these partial decay rates 
through the corresponding reduced decay rates (|E.ip . For that we note that for an A throat KK mode, we have 
jTV'^pZA ~ 1 and \N^\'^zs ^ Prszs / za, where the tunnelling probability Prs ^ m'^/k'^. We therefore have 



which, together with ljE.6p . implies 



Ta 



za 



zs 



SX' 



2S^ 



ZA ITl 

Repeating the same argument for an SM throat KK mode leads to an analogous equation. 



^oPrs 



ZA 

Zs 



SX' 



2S_ 

m 



s 



Here 6X^ is given by an expression analogous to (|E.5p . namely 

5S 



5X- 



zs 



(E.7) 



(E.: 



(E.9) 



(E.IO) 



We now have two equations, (|E.8(I and IjE.Qp . for the two unknown sinks. (Note that while we have written 
the equations up to order one factors, one could write them exactly, though in a much more cluttered and less 
intuitive form.) The solution can be written as 



m 



1 



Pns^ 
Zs 



2 1 , m{zs-ZA) 



1 



Pns^ 

ZA 



PRS^)-rl(l-Pns^ 

zs J \ ZA 



(E.ll) 



Clearly, for a large enough barrier height k/m (such that the factors Pb.szs/za and m{zs — za) / i/^k"^ zsza) 
are both small), we find ~ F^m/2, - F^m/2 as expected; the presence of the finite barrier induces 
corrections suppressed by various powers of the barrier factor. 
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